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X-ray Hardness Evolution in GRB Afterglows and Flares: 
Late Time GRB Activity Without Nh Variations 

Nathaniel R. Butler 1 ' 2 and Daniel Kocevski 2 
ABSTRACT 

We show that the X-ray and 7-ray spectra of Swift GRBs and their after- 
glows are consistent with the emission characteristic of an expanding, relativistic 
fireball. The classical afterglow due to the impact of the fireball on the external 
medium is often not observed until one to several hours after the GRB. Focusing 
on GRBs 061121, 060614, and 060124, but generalizing to the full (>50 Msec 
XRT exposure) Swift sample up to and including GRB 061210, we show that 
the early emission in >90% of early afterglows has a characteristic vF v spectral 
energy E pea k which likely evolves from the 7-rays through the soft X-ray band 
on timescales of 10 2 — 10 4 s after the GRB. The observed spectra are strongly 
curved when plotted with logarithmic axes and have often been incorrectly fit- 
ted in other studies with a time-varying soft X-ray absorption. The spectral 
evolution inferred from fitting instead models used to fit GRBs demonstrates a 
common evolution — a powerlaw hardness intensity correlation and hard to soft 
evolution — for GRBs and the early X-ray afterglows and X-ray flares. Combined 
with studies of short timescale variability, our findings indicate a central engine 
active for longer than previously suspected. The GRB spectra are observed to 
become very soft at late times due to an intrinsic spectral evolution and due 
to the surprising faintness of some afterglows. We discuss models for the early 
X-ray emission. 
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Introduction 



The Swift satellite (IGehrels et al.l 12004 ) and its X-ray telescope (IBurrows et al.ll2005bl ) 
have opened a new window into the early lives of 7-ray Bursts (GRBs) and their af- 
terglows. We see a complex array of behaviors, many of which appear to directly con- 

the well tested 



flict (e.g., lOTjrien et al.ll2006l : IPanaitescu et al.ll2006 



Willin gale et al 



pear t 
l2006h 



1994; Sari fc Piran 



internal-external shock GRB and afterglow model (Rees fc Meszarosl ^ 
19971 : ISari. Piran. fc NaravarJ Il998l : Iwiiers fc Galamal Il999h . In this "fireball" model, the 
GRB is produced via collisions of shells in a relativistic outflow, and an afterglow arises later 
as the ejecta sweep up and heat the surrounding medium. The Swift afterglows exhibit dra- 
matic flaring, rapid ly decaying prompt emission tails, and typically a broad plateau phase 

'ly afterglow observations prior to Swift (e.g., 
10s duration burst rapidly gone and replaced 



until t « 10 4 s (e.g. 


Nousek et al. 


2006) 


Frontera et al. 


2000 


) suggested instead 



by the fading afterglow emission. How these observations are to be reconciled and what 
mechanisms produce the early afterglow emission are key open questions. 

Particularly intriguing, several recent stu dies fit the Swift X-ray Telescope (XRT) data 
and in fer a time variable soft X-ray absorption ( Starling et al. 2005 ; Rol et a02006 ; Campana et al 
2006a ). This would imply that the early afterglow is stripping electrons from a dense shell 
of light-element-rich material located R <^ 1 pc from the GRB, which was not already fully 
ionized by the GRB. It is difficult to detect su ch an effect because of the strong spectral 



evolution common in the early afterglows (e.g., iVaughan et al.ll2006l ; iButlerl l2007al . "Paper 
I" ). A changing column density Nh cannot easily be separated from intrinsic afterglow spec- 
tral evolution, given the narrow XRT bandpass. If the early X-ray spectra exhibit log-log 
curvature l ike that of GRBs, which have vF v spe ctral turnovers at a characteristic energy 



-E-peak (e.g., iPreece et al.ll2000l ; iKaneko et al.ll2006l ). then evolution in the curvature could be 



mistaken for variations in Nh- 

As we discuss below, plots of early XRT spectra do show strong log-log curvature and 
an inferred -E pea k which typically passes in time through the X-ray band. This produces 
a changing X-ray hardness, which we observe to correlate with the flux. A close anal- 
ogy can be found in the spec tral evolution of GR Bs observed by the Burst and Transient 
Source Experiment (BATSE; iFishman et al.lll999l ). A characteristic feature of these spec- 
tra and ligh t curves is a hard-to-soft evolution in time and a powerlaw hardness-intensity 
Golenetskii et alill983l: iKargatis et~al~1ll995l : iNorris et allll996l: iFenimore et al.l 



corre 



1995 



at ion 



Fenimore. Madras, fc Nayakshin 19961 ). The recent refined study of Borgonovo fc Ryde 



20011 ) measures a powerlaw relation between the characteristic energy -Ep ea k and the bolo- 
metric flux Fboi valid for >57% of GRB pulses, -Ep ea k oc F®^ ±0 ' 2 . We observe a consistent 
correlation in the soft, early-time XRT data. 
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In Paper I, we present evidence for this outlier population of extremely soft afterglows 
in the first year of Swift XRT afterglow data. Although they were identified via an auto- 
mated search for spectral lines, the spect ra are also fitted we ll by m odels containing multiple 
continuum components. Below and in iButler &: Kocevskil (120071 ). we explore further the 
phenomenology associated with this soft emission. We demonstrate that GRB-like behav- 
ior is present in the first t <^ 1 hour of >90% of the afterglows and is especially promi- 
nent during the flaring. In two cases, thanks to Burst Alert Telescope (BAT) triggers on 
bright precursors, X-ray emission coincident in time with the classical GRB is detected and 
can be shown to have quite similar properties to the highly time-variable emission at later 
times. This is strong evidence; — to be combined with the short timescale variabiliity stud- 
ies (e.g., 



Burrows et al.ll2005al ; iFalcone et al.ll2006l ; iRomano et al.ll2006al ; iPagani et al.l 12006 



Kocevski fc Butler 
engine. 



20071 ) — tying the flare and early afterglow emission to the GRB central 



2. Data Reduction 

Our automated pipeline at U. C. Berkeley downloads the Swift data in near real time 
from the Swift Archivqj and quicklook site. We use the calibration files from the 2006-04- 
27 BAT and XRT database release. The additional automated processing described below 
is done uniformly for all afterglows via custom IDL scripts. The final data products are 
available for general consumption^. 

The XRT suffers from a significant number of bad or unstable pixels and columns. Two 
central detector columns were lost due to a micro-meteorite strike^. For the early afterglows 
(t <; 10 3 s), when the satellite initially points the XRT at the source without the precise 
localization information needed to offset from the bad columns, a large and time-dependent 
fraction of the flux can be lost. In order to produce accurate light curves and properly 
normalized spectra, it is necessary to accurately determine the position centroid and to 
precisely track the loss of source and background flux due to the bad detector elements on 
short (~ few second) timescales. 



: ftp: / /legacy.gsfc. nasa.gov/swift /data 
2 http: / / astro, berkeley.edu/~nat/swift 

3 http://swift.gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/SWIFT-XRT-CALDB-01_v5.pdf 
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2.1. Photon Counting (PC) Mode Light Curves 

We begin by projecting the data in the 0.5-8.0 keV band from each PC mode followup 
observation onto a tangent plane centered at the source position quoted by the XRT Team. 
In raw coordinates, we reject all pixels with more than six counts and also containing more 
signal than contained in the surrounding 8 pixels summed. Using the aspect solution file 
(*sat*.fits), we determine the satellite pointing for each detection frame. We then map the 
bad pixels in raw detector coordinates determined by xrtpipeline and by our algorithm 
onto the sky on a frame-by-frame basis. This is used to generate exposure maps for the full 
observation and function of time. 



Using the full exposure map, we determine the afterglow position centroid (see, [Butler 



2007bl ) to fix the source extraction region. We consider a 16 pixel radius source extraction 



region, surrounding by a n annular bac kground extraction region of outer radius 64 pixels. 



Running wavdetect (see, iButlerl l2007bl ) . we then determine the positions of field sources in 
the image. We mask out the regions corresponding to the field sources from the source 
and background extraction regions. Also, using the Point Spread Function (PSF) model 
(swxpsf 20010101v003 . f its) at 1.3 keV, we determine the level of residual field source 
contamination in the source extraction region (typically negligible) for later subtraction. 

Initially ignoring pileup, we extract the source and background counts for each good 
time interval of data acquisition. The fraction of lost signal and the scale factor relating 
the background in the source and background extraction regions is determined for each 
extraction using the time-dependent exposure map. assuming these exposure corrections for 
the entirety of each time interval, we subdivide the counts in each interval so that a fixed 
signal-to-noise of 3 is achieved. 

In order to check and to account for pileup, we perform a coarse Bayesian blocking 



(jScargkj Il998l ). with a strong prior weight against adding a new segment (e~ 50 ). Using 
the maximum observed count rate in each segment thus determined, we find the minimum 
aperture necessary to reduce the source signal to levels where pileup is negligible. The coarse 
blocking results in a small number of regions (typically 2-3) of differing inner extraction 
radi us for an afterglow . We assume pileup is important for count rates > 0.5 cps (see 



also, iNousek et all hood ) . The light curves are verified to transition smoothly across regions 



of different inner extraction region radius. Using the time intervals and pileup corrected 
apertures thus determined, we rebin the data to a signal-to-noise of 3 and recalculate the 
exposure correction for each time interval. The final time regions and exposure corrections 
define our temporal extractions regions for the extraction of light curves in different energy 
bands and for extraction of spectra below. 
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2.2. Windowed Timing (WT) Mode Light Curves 

Our reduction of the WT mode data closely parallels our PC mode reduction, except 
that it is more natural to extract the WT mode data in raw detector coordinates than in sky 
coordinates as done above for the PC mode data. This is due to the readout mode; detector 
pixels are summed in RAW-Y and the resulting data are in column (RAW-X) format. 

Summing the data from each WT mode followup, we reject any RAW-X columns con- 
taining a > lOcr count rate relative to the background, after first ignoring pixels in the 16 
pixel source extraction region. We also reject any RAW-X columns containing 100 times 
more signal than the highest neighboring column (or > 100 if the neighbors contain no sig- 
nal). Using the sky image determined from the PC mode data and the satellite aspect, we 
project the background onto the RAW-X axis and form a background mask for the 64 pixel 
outer radius and 16 pixel inner radius extraction region. We do not allow masking of the 
pixels within the central 16 pixel source regions. If the source is bright (> 10 3 cps), we re- 
center the source and background apertures. Small aspect shifts ~1 pixel are not uncommon 
between the PC and WT mode data and must be accounted for. 

We determine the exposure corrections as for the PC mode data, but also adjusting 
the PSF model for the WT mode summing of RAW-Y pixels. We note that our exposure 
corrections account for source signal contained in the background region. We determine 
a pileup correction as above, but with a limiting source count rate of 150 cps (see also, 



Nousek et all 120061 ). 



2.3. PC and WT Mode Spectra 

Spectral response files are generated using the xrtmkarf task for each time interval of 
interest. Our invocation of the task ignores the exposure maps calculated above, determin- 
ing the energy dependence of of the source extraction assuming only the inner and outer 
source and background regions. We then adjust the normalization of the resulting Ancillary 
Response File (ARF) to account for the actual loss in flux (0.5-8.0 keV) on a pixel by pixel 
basis using the divided time intervals and associated exposure corrections determined above. 
The spectra are fit in ISIsfl For each spectral bin, we require a S/N of 3.5. We define S/N 
as the background-subtracted number of counts divided by the square root of the sum of 
the signal counts and the variance in the background. As done in Paper I, we restrict our 
attention to time-resolved spectra containing 500 or more counts or to spectra formed by 



4 http: / /space.mit.edu/CXC /ISIS / 
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grouping two or more of the 500 counts spectra. 

We fit the PC and WT mode data over the 0.3-10.0 keV range, also accounting for the 
systematic calibration uncertainties ~ 39c(E In WT mode, we allow the detector gain to vary 
by ±80 eV@. 

2.4. BAT Light Curves and Spectra 

We establish the energy scale and mask weighting for the BAT data by running the 
bateconvert and batmaskwtevt tasks. The mask- weighting removes flux from background 
sources. Spectra and light curves are extracted with the batbinevt task, and response 
matrices are produced by running batdrmgen. We apply the systematic error corrections 
to the low-energy BAT spectral data as suggested by the BAT Digest website^, and fit 
the data using ISIS. The spectral normalizations are corrected for satellite slews using the 
batupdatephakw task. For GRB 060124 below, BAT spectral fits are performed on the mask 
tagged light curve data in four channels, assuming the on-axis response and also accounting 
for the systematic error. 

3. The Joint BAT+XRT Spectra of Three Events 

There are two bright events in the XRT sample which overlap in time entirely with what 
would commonly be thought of as the prompt phase of GRB emission. The observations by 
the XRT were made possible by a bright precursor just minutes prior to each GRB observed 
in the BAT, on which the BAT triggered. We therefore have both BAT and XRT data for 
each event, GRB 060124 and GRB 061121. We also discuss the bright event GRB 060614, 
which has excellent XRT coverage due to an early, rapid spacecraft slew. 

Figure [1] displays spectral fits to a selected set of time-resolved intervals in each events. 
The best-fit model parameters are given in Table 1. The time evolution of these parameters 
are presented and discussed in detail in the next three subsections. 



5 http: / / swift.gsfc.nasa.gov/docs /heasarc/caldb /swift / docs /xrt / spie05_romano.pdf 
6 http: / /swift. gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/xrt_bias.pdf 
' http: / / swift.gsfc.nasa.gov/docs/swift /analysis /bat _digest.html 
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3.1. GRB 060124 



Swift-BAT triggered and located the precursor to GRB 060124, allow ing the XRT to 
slew and begin s imultaneous observatio ns 106s later ([Holland et al.l 120061 ) . This event is 
also discussed in iRomano et al.l (l2006bl ). The 0.3-10.0 keV light curve is plotted in Figure 
[2J There are two prominent peaks. As shown in the background (lighter two shades of 
gray), the time profile in the soft XRT channel (0.3-1.3 keV) is broader than that in the 
hard channel (1.3-10.0 keV). The BAT light curve shows even narrower time structure and 
resolves the broad first XRT peak into at least 3 sub-peaks. The light curve after the flare 
(t > 10 4 s) and extending to 22 days is well fit by a powerlaw t -1 - 32 ^- 01 (x 2 /v = 535.2/465). 

We group the XRT data into ^ 500 counts spectra and fit powerlaws (Figure [3], left). 
Each fit is statistically acceptable, with a reduced x 2 °f order unity. The photon index T 
is observed to decrease in time, although with modulation in time that correlates with the 
X-ray flux and with Nh (see explanation in Section l5"TTj) . At late times (t > 10 4 s), the Nh 
values asymptote to the blue, dashed curve (N H = 2.3 ± 0.2 x 10 21 cm -2 ) plotted in the 
figure. 

To study the time varying log-log curvature, we jointly fit the BAT and XRT data using 
the iBand et al.l (119931 ) model. Here, we choose extraction regions which allow for a BAT 
signal-to-noise of 20 or higher. We also fix the column density Nh to the late time valu e. The 



model fit is actually a progression of fits of nested models (e.g.. iProtassov et al.ll2002l ). from 
the simplest powerlaw model to a powerlaw times exponential model, to the smoothly broken 
powerlaw Band model. Each more complex model has one additional degree of freedom. We 
accept or reject the more complex model at each stage by requiring A% 2 > 2.706 (i.e., 90% 
confidence). If the data are acceptably fit by only the powerlaw model, we quote a limit 
on i^ak using eit her the exponential ti mes powerlaw model (for T < 2) or the constrained 
Band formalism (jSakamoto et al.l 120031 . ; for T > 2). In order that -E pea k correspond to a 
peak in the vF v spectrum, we require the low energy index a > —2 and the high energy 
index (3 < — 2. After finding tha t the fits were consi stent with a < 0, as also found for 
BATSE GRBs JPreece et al.lboool : iKaneko et aDbood ). we included this as a constraint to 
derive the tightest error bounds on the other model parameters. 

As shown in Figure [3] (right), the data are better fit (> 90% confidence) with the Band 
model in most of the time regions. The peak energy rises and declines with each of the 
four prominent light curve pulses. For each pulse, we present powerlaw fits to the E veak 
decli nes. The rises are not well me asured, as is also typically the case for BATSE bursts 



[e.g.. iKocevski. Ryde. fc Liang! 120031 ). Prior to t ^ 800s, the observed spectrum corresponds 
mostly to the low energy portion of the Band model spectrum, except episodically at the 
flare troughs, where E^^ enters the X-ray band. These times regions are also those of 
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highest Nh in Figure [3] (left). The third pulse decline exhibits a strong evolution in both a 
and -Epeak- After t ~ 800s the observed spectrum corresponds to the high energy portion of 
the model spectrum, and -E pea k has transited the X-ray band. Figure [T] (middle) plots the 
vF v spectrum at 3 time epochs. 



Motivated by the watershed event GRB 060218 (Paper I; ICampana et al.l l2006al ). we 
also attempt to fit the X-ray curvature using a powerlaw plus blackbody model. The fits to 
the X-ray data alone are provocative and show a smooth temperature decline after each of 
the two major pulses. However, the fits are statistically unacceptable when we also attempt 
to account for the BAT data. This is also true for the GRB 061121 spectra discussed in the 
next sub-section. This should be taken as a caveat also to the powerlaw plus blackbody fits 
presented for the XRT data in Paper I, where the derived blackbody temperature variation 
may imply instead to -E pea k variations. We note, however, that the X-ray spectra of the 
unusual GRB 060218 burst and afterglow are better fit by a blackbody plus powerlaw than 
by a Band model (Paper I). We do not consider the possibility of two powerlaws and a 
blackbody for the bursts discussed here. 



3.2. GRB 061121 



Swift-B AT triggered on and began observing the pre cursor to GRB 06 1121 55s prior to 
the XRT slew toward and onset of the mai n GRB event (IPaee et al.ll2006l ). The Swift team 
designated this event a "Burst of Interest" JGehrels et alJbood ) due to the rare simultaneous 
detections in the BAT and XRT bands and at longer wavelengths. As shown in Figure HJ the 
7-ray and X-ray light curves show multiple peaks, with most of the prominent time structure 
apparent in only the 7-ray band. 

In Figures [3] (top), we show the results of powerlaw and Band model fits to the 061121 
data. The data are not of as high signal-to-noise in the X-ray band as the 060124 data, 
however, many of the same trends are apparent. There is a hard to soft evolution apparent 
in the powerlaw photon index and a correlation between the photon index and Nh- The 
Band model photon index goes from the low energy side to the high energy once E pea k has 
crossed the X-ray band. E pea k also appears to rise and fall with flaring prior to 80s. The 
vF v spectrum is plotted at two epoch in Figured] (top). For the Band fits, we use the late 
time (t > 10 4 s) N H = 2.5 ± 0.3 x 10 21 cm" 2 . 

XMM da ta for this event be ginning; after t « 6hrs show consistent powerlaw fits with our 



late-time fits (IButler et al.ll2007l ). In particular, Nh = l>71_ao2 x 10 cm , consistent with 



our late-time Nh at the 2a level and well below the mean early-time value. XMM and XRT 
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data generally agree well with respect to the late-time Nh determinations (e.g. jMoretti et al. 



2006|). 



3.3. GRB 060614 



The GRB 060614 (IParsons et al.ll2006l ) afterglow fades rapidly as a powerlaw from the 
prompt emission, with no flaring (Figure [5]). There is excellent BAT+XRT coverage during 
the prompt tail emission lasting to t ~ 150s. We observe weak Nh-T correlated modulations 
during the rapidly fading tail, which would imply an Nh that decreases in time, reaching 
the value marked by the dashed line in Figure [3] (left) by t = 10 4 s. 

However, the Band model fits show an -Ep ea k which passes through the X-ray band 
without requiring a varying Nh- Extrapolating backward through the prompt emission, the 
best fit decay also fits two -Epeak limits derived for the BAT only prompt emission. Figure 
[T] (bottom) plots the vF v spectrum at two epochs. Expressed in terms of the flux F X rt as 
measured by the XRT rate, E peSLk oc F-° T 72±a03 . 

The low-energy photon index a also appears to evolve in time after the main GRB 
emission. 



3.4. Hardness Plots for GRBs 061121, 060124 and 060614 

It will be useful below to see how the spectral evolution in the early X-ray light curves 
of GRBs 061121, 060124 and 060614 impacts the X-ray hardness ratio. We define this is the 
ratio of counts in the 1.3-10.0 keV band to the counts in the 0.1-1.3 keV band. The average 
hardness ratio {HR) for most afterglows is 1. 

Figure[6]in 9 panels shows the hardness and rate time profiles for GRBs 061121, 060124, 
060614. The middle panels (looking top to bottom) sho w the X-ray light curve fit using an 



extension of the Bayesian blocks algorithm (IScarglelll998l ) to piecewise logarithmic data. The 



rate and hardness data are fit jointly, allowing the minimum number of powerlaw segments 
such that x 2 l v ~ 1- The fits to the rate and hardness are plotted in the top and middle 
panels, indexed according to time. The hardness tracks the flux and moves along roughly 
parallel tracks. In the bottom panels, the flux in both XRT bands (top panel) and the 
hardness (bottom panel) are plotted for each powerlaw segment. During the decline phase of 
each pulse, the hardness scales as the square-root of the rate for GRBs 061121 and 060124. 
For the GRB 060614, the hardness and flux track as found above for E pea ^ and flux. 
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Each pulse in GRB 060124 peaks at roughly the same time, independent of energy band. 
There is, however, a hardness rise during the flux rise because the hard band increases more 
rapidly. There is also a modest overall hard to soft trend throughout the light curve. 

The hardness plot does not capture the strong spectral variations between 500 and 600s 
in GRB 060124, which are apparent from the broad band fits (Figure [3] middle) and occur 
mostly for E peak above the XRT bandpass. The time dependences of E peak during this region 
and later are given in the figure. The -E pea k dependence can also be given in terms of the flux 
F, in order to sidestep the problem of unknown start time. For all but the last flare, where we 
use the XRT count rate, we use the BAT 15-350 keV count rate for the flux. For pulses 1-4, 
we find E vcak oc F^ ±17 , i^ T 8±a5 , F^ ±0 -\ F^ ±0 - 2 . In the bottom right left panels of 
Figure El we show that the hardness can be described by the square root of the observed flux, 



as is common for GRBs at higher energies observed with BATSE (e.g.. iBorgonovo fc Ryde 



200ll : iRvde fc Petrosianl 12002 : iKocevski. Rvde. fc LianeJ 120031 : iRvdel l2005l . Section E> . 



For GRB 061121, the hardness plots show an initial hardening followed by a decrease in 
the hardness which scales well with the square root of the X-ray rate. There may be broad 
pulses on top of the decline, although these have only a minor impact on the hardness. 
GRB 060614 appears to mostly to exhibit a secular decline in both flux and hardness, 
corresponding to the fading tail of the prompt emission. 

For each GRB, the hardness plot capture the Speak evolution in general terms. Both 
(HR and -Epeak) decrease during rate declines at a similar power of the rate. It is apparently 
not possible to cleanly if at all separate evolution of a from evolution of -Epeak, given the 
hardness alone. From Figures [3] (top right) and [3] (middle right) and also from time-resolved 
spectral studies of many GRBs (Section [5]), these parameters tend to evolve simultaneously. 



4. Example Spectra for 4 Other Early Afterglows 

Most early X-ray afterglows have a low signal-to-noise or no coincident detection by the 
BAT. It is possible to derive -E pea k values or limits for these early on, given the BAT data. 
Late time Epeak from the X-ray data typically show values in or passing through the XRT 
band after one to se veral minutes. The spectral evolution from one such event, GRB 060714 



( iKrimm et al.ll2006l ). is shown in Figure 



The hardness plot (Figure [7] middle) allow for a finer time sampling of the spectra 
evolution. The hardness (likely also E pea k) rises and declines with the flux along the same 
track in the hardness-rate plane as two flares take place. The column density (not plotted) 
is a factor two larger in the time interval 140-1 70s than outside that interval, indicating an 
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Epeak passage. There are a handful of examples with higher signal-to-noise XRT observations. 



The GRB 060526 (jCampana et al.ll2006bl ) afterglow exhibits time-correlated T-Nh vari- 
ations and a corresponding rapid then smooth decline of -E pea k through the XRT band (Figure 
[8] left). The initial GRB pulse (t < 9.4s) is well fit by a simple powerlaw (a = —1.6 ± 0.2, 
X 2 /v = 16.83/16), and we derive £ pcak > 80 keV (90% confidence). The flare at t ~ 250s 
is detected by the BAT as well, and we use the BAT data to obtain the best Band model 
fits. The Band model photon indices are poorly measured. The composite flare and decline 
is shown in Figure [9J The hardness evolves similarly to the best-fit -E peak values. 



The very bright afterglow to GRB 060729 (IGrupe et al.ll2006al ) continues to be detected 
4.5 months after the GRB. The GRB is over and done with by t ~ 130s in the BAT. We 
find -Epeak > 50 keV (90% confidence). After t > 100s in the XRT, there is a rapid decline, 
interrupted by a flare or rise at 160s (Figure [9]). Time-correlated N^-Y variations and an 
-Speak passage through the X-ray band are similar to those discussed above (Figure [8]). We 
observe that -Epeak declines with the X-ray rate as 0,1 both before and after the mild 

flare at t ~ 180s. There is also a possibly significant decline in j3 with time. 

The hardness declines by an order of magnitude, reaching a minimum at t ~ 250s, 
and then increases to the late time (t > 10 3 s) value. Note that no clear coincident change 
is present in the rate plot. The hardness plot demonstrates that the late time emission is 
spectrally different from the early emission and that its onset occurs at t ~ 250s. 



Modest but clear Nh~T variations are seen for GRB 060904B (IGrupe et al.ll2006bl ). The 
prompt emission [t < 8.3s) has Epeak = I25I30 5 keV. -E pea k transits the X-ray band nicely 
(Figure [8]). The hardness evolution shows the usual time dependence in the declining tail of 
the flare (Figure [9]). E peak decays versus the rate as F^^ ±0 2 . 



The emission for GRB 060929 flMarkwardt et all 120061 ) at t < 13s exhibited £ peak > 75 
keV. The X-ray flare peaking at t ~ 550s is weakly detected by the BAT. In the XRT, there 
is a clear softening trend (Figure [9]), likely Nh-T variations, and an -E peak declining through 
the X-ray range (Figure [S]). -E peak drops with the X-ray rate as ^xrt ■ The hardness 
reaches a minimum at t = 630 ± 10s. 



5. Discussion 

5.1. Global Sample Properties 

In terms of the spectral evolution properties, we see no apparent difference between the 
fading tales of flare-like X-ray emission and the rapid X-ray declines often observed to trail 
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flaring in the BAT (e .g. , iTagliaferri et al.ll2005l ; iBarthelmy et al.ll2005l ; ICusumano et al.ll2006 



Vaughan et al.ll2006l ). Indeed, based solely on timing properties, many of the rapid declines 



rapid declines are thought to be the fading tail of the prompt emission 


Panaitescu et al. 


2006; 


Yamazaki et al. 2006; Lazzati & Begelman 2006; Zhang et al. 


2006 


), and the X-ray 


flares are thought to be due to later central engine activity ( 


Zhang et al. 


2006: Ioka et al. 



20051 ; iFan fc Weill2005l ). We observe a clear distinctions between the spectra measured before 
the light curve plateau and after the start of the plateau; only t he lat e spectra exhibit a tight 
clustering with r « 2 (Figure [TO] Paper I; iButler &; Kocevskil (120071 )). 



Figure [TT] shows what we expect to measure from powerlaw fits to a time-evolving Band 
model spectrum. As -E pea k enters the X-ray band, the spectral curvature as would be seen 
on a plot with logarithmic axes increases and the inferred X-ray column density increases 
linearly with an increasing inferred photon index T. This occurs despite the fact that only 
-E'peak changes in the simulation. Figure [12] suggests that the effect is common in the XRT 
data (Section 15.21) . 

Figure d3] (left) shows that the flares (Table 2) and rapid X-ray declines exhibit signif- 
icant hardness-intensity and hardness-fluence correlations which match closely the correla- 
tions observed for GRBs (Section 15.31 below) . 



; is common 


:o observe finer time structure at higher energies as compared 


(Norris et al. 


1996: 


Fenimore et al. 


1995: 


Fenimore. Madras. & Navakshin 



19961 ). Pulses tend to be narrower, fade more rapidly, and evolve stronger spectrally at high 
energies. Consistent with this, the X-ray flares (and also the rapid declines) appear longer 
(8 ± 1% on average, Figure [Til left) and with smoother time structure (e.g., Figure [2]) at 
softer energies. This can be understood as the effect of E peak evol ving into the X-ray band, 



which allows the X-ray emission to be observed for longer (e.g., iKocevski. Ryde. fc Liang 



20031 . and Section [53]). Although it is difficult to see by a eye, there is also evidence for a 
25 ± 5% increase in the flare rise time with decreasing X-ray energy band (Figure [Til right) . 
This is close to the expected pulse broadening fraction from an extrapolation of the GRB 
behaviour, 1 — (1.3/0.5) -0 ' 4 ~ 30%, where 0.5 and 1.3 keV are used as approximate lower 
bandpass energies. Given the possibility that resolved 7-ray flares are blurred together in the 
X-ray band (e.g. Figure [2]), however, it is not clear how meaningful this apparent consistency 



is. 
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5.2. The Physical X-ray Column Density Does Not Vary 



The time-resolved XRT afterglows are well fit by absorbed powerlaws at all epochs (see 
also Paper I). Prior to a characteristic hardness variation turn-off time T h ~ 10 2 — 10 4 s, 
which we discuss for a large sample of bursts in (IButler fc Kocevskil 120071 ). there is strong 
evolution in both the best-fit photon indices and the best fit column densities Nh- After this 
time, the quantities typically do not vary. To fit more complicated models to the early time 
afterglows, we have found it necessary to jointly fit the BAT and XRT data (when possible) 
and to tie the column density to the value measured at late time. The late time value is 
typically not the Galactic value. 

Band model fits are able to account for both the BAT and XRT emission without a time 



variable column density (see, also. iFalcone et al.ll2006l ). The ubiquitous hardness evolution 
appears to be best understood in terms of an evolving E peab ^ as we discuss in detail below. 

Several studies have claimed recently a decreasing N H b ased on fits to the XRT data 
Jstarling et all 120051 : IRoI et al.l I2OO6I : bampana et~alTl2006cl . GRBs 050730, 050716, and 
050904, respectively). Each study presents a coarsely time- resolved set of spectral fits, which 
demonstrate a higher Nh at early times. This is an artificial feature that we observe in fits 
to most Swift early afterglows. It is especially clear in the brightest afterglows, which often 
sample the declining tail of the prompt emission. For each of the 3 bursts with claimed Nh 
variations (e.g., Figure [T5|) . a fine times-scale spectral analysis reveals an Nh which both 
increases and decreases in time (following T and the flux). Observed drops in Nh ^ 10 21 
cm -2 (or ^ 10 22 cm -2 in the rest-frame) on timescales of 10 — 100s are challenging enough, 
but drops and increases and drops again on these timescales are unphysical. 

We strongly caution against taking the early Nh values at face value. Measurements of 
Nh at t <z 10 4 s will be artifically high. Also, although we cannot rule variations out in all 
cases, they are not required by the data and they are also not the simplest interpretation of 
the data. Firm measurements of Nh variability will require finely time-sampled broadband 
data (e.g., Ultra-violet, X-ray, and 7-ray data) to disentangle the effects of the evolving Band 
model spectrum from the soft X-ray photoelectric absorption. 

For those fitting XRT spectra, we recommend measuring Nh at late times (t ^ 10 4 s) or 
performing joint fits at different time intervals with a single Nh parameter shared between 
multiple spectra. Jointly fit with BAT when possible. Fine time resolution is essential 
when testing varia ble Nh', it is not sufficient to fit expo nential times powerlaw models or 
Band models (e.g.. lRol et al.ll2006l ; ICampana et al.ll2006a ) with coarse time resolution. The 
hardness ratio can be utilized to diagnose cases where inferred Nh values are likely to vary 
artificially. 
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5.3. Is the Early XRT Emission the Same as Prompt GRB Emission? 



We have shown for seven events that the early X-ray spectra requ i re a fit model which 



also has been shown to reliably fit a 



Frontera et al. 2000; Sakamoto et al 



1 GRBs (e.g., IPreece et al.ll2000l : iKaneko et all 12006 



200 3]). The need for such a model is also clear from 



hardness variations (see, also. lButler &: Kocevskill2007l ) and time correlated iV#-r variations 
observed for even low signal-to-noise afterglows, which demonstrate a characteristic increase 
in spectral curvature in the XRT band. 

In cases where E peilk is well measured, or using HR when -E pea k is poorly measured, we 
observe a hard-to-soft evolution and a strong hardness-intensity correlation, also commonly 
seen in GRB pulses. Our correlation can be described as a hardness which tracks the flux 
to a power 0.43 ± 0.07. From the Band model fits, our best fit E pea ^-F relation is 0.7 ± 0.2 
(Table 3). A closely consistent powerlaw relation exists for most GRB pulse s, also with a 
large scatter in observed values (jGolenetskii et al.lll983l ; iKargatis et al.lll995l). T he scatter 
is apparently minimized for bolometric measures of flux (jBorgonovo &: Ryddl200ll ). yielding 



-E'peak oc F° o j ±0 ' 2 . The fact that we have observed a consistent relation can be turned around 
to imply GRB-like emission with -Z? pe ak ~ IkeV, typically. That Swift observes bright X- 
ray flares appears to be a consequence of this and also due to the surprising fact that the 
afterglow is faint at these times. It is interesting to speculate that there may be bright 
optical flares due to internal shocks at times of several hours after some GRBs with faint 
afterglows. 

The typical E ppa k values for the XRT are two orders of m agnitude below the mode of the 
BATSE distribution (IPreece et al.ll2000l ; IKaneko et al.ll2006l ). As we discuss below, some of 
the soft E peak values may be due to viewing effects of delayed emission with an intrinsically 
higher E pe3 k. However, the soft flare emission implies intrinsic spectral evolution or soft late 
central engine activity which would extend the BATSE E pea ^ distribution. Our derived values 
for a are poorly constrained, but likely consistent with the BATSE distribution. Finally, it 
is remarkable that very soft emission is observed in a few cases, extending the distribution 
in (3 to very low values < —6 (Figures M and [TDJ GRBs 050714B and 050822 discussed in 
Paper I). 



5.4. Interpretation of the Spectral Variations 



Although intrinsic spectral evolution is likely also present, most of the softening trend 
and h ardness-intensity correlation in GRB pulses is attributed to the so called "curvature ef- 
fect" ( Fenimore. Madras, fc Nayakshin 1996 ; Sari &; Piran 1997 ; Norris|[2002 ; Ryde fc Petrosian 
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20021 : iKocevski. Rvde. fc LianeJl2003l : bin et alJfiooi : bin fc Lull2005l : Ishen. Song, fc Lill2005h . 



(Nousek et al. 


2006: 


Zhang et al. 


2006: 


Panaitescu 


2007) 



Derivations from first 



princi p les of the curvature effe ct on the observed spectra can be found in lGranot. Piran. fc Sari 
f)l999h : IWoods fc Loebl fjl999h . 



If we imagine a spherical emitting shell at radius R that emits as a delta function 
at t , the spectral flux F E scales with the Doppler factor 8 as F E oc F E [E5}/5 2 . Here, 
5 = 7(1 — /3 c cos(#))/(l — /3 C ) pa 1 + 7 2 # 2 , where 9 is the viewing angle to emitting material 
off the line of site. The photons from larger angles will be delayed in time t — t Q = (1 + 
z)(6 - 1)(1 - Pc)R/(cP c ) » (1 + z)6 2 R/(2c) oc 5. 

For a powerlaw spectrum F E oc E 1 ^^, the observed fl ux declines in time as a pow erlaw 
(t — t ) - ' 6 , with |6| = 1 + |a| and no hardness evolution ( Kumar &: Panaitescul l2000l ) . For 
a Band spectrum, we see either the low energy index a or the high energy index f3 or some 
average of the two, depending on the location of -E pea k with respect to the bandpass. -E^k 
will decline as (t — £ ) _1 - When -E pea k is in the band, the vF v turnover implies — a e fr pa 1 — 2, 
and we expect to see a powerlaw hardness-intensity correlation -Ep ea k oc _p a3 -°- 5 . Larger 
values of the index are favored observationally, because they correspond to a higher flux. 
We will observe the hardness (which our simulations show to scale linearly with -Epeak for 
a range of Band model parameters) to approximately linearly correlate with the fluence. 
Departures from this expected behavior will occur for emitting shells of different shape, for 
an inhomogeneous emitting surface, for non- instantaneous emission, or if intrinsic spectral 
evolution domin ates. Also, the m easured flux decay in time is a strong function of the 



assumed t (e.g.. lLiang et al.ll2006l ). 



Our best fit HR-F relation index (Figure [T3"j) and our average E^^-F relation rela- 
tion index (Section I5.3t Table 3) are consistent with those expected in this simple picture. 
Spectral variations are n ot inc onsistent with the curvature effect, as recently suggested by 
Zhang. Liang. &: Zhang! (120061 ). Rather, they facilitate a higher order test of the curvature 
effect, and allow us to confirm the curvature effect in way that shows the X-ray phenomenol- 
ogy to closely parallel the 7-ray phenomenology. Moreover, the scatter in our HR-F relation 



(Figure [T31 left) is less than that found for time- index-energy index relations (INousek et al. 
20061 : iPanaitescul 120071 ). which assume powerlaw X-ray spectra. 



The mean time index for the Speak decays in Table 3 is — 1.4±0.6, consistent with unity. 
This indicates that our choice to associa te t„ with the star t of the flare or pulse is roughly 
correct, in agreement with the findings of iLiang et al.l (120061 ). Although we see evidence that 
later flares often have lower -E pea k in the same event with multiple flares (e.g., Figure [3]), we 
do not see a correlation between t in Table 3 and -E pea k just after that time. 
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We have observed two cases of a evolution (Figures[3]right) which accompanies the E vea ^ 
evolution. Due to the proximity of Speak to the bottom of the XRT pass band and also due 
to the possibility of a modestly incorrectly measured N H , these cases should be interpreted 
cautiously This evolution, or that observed for f3 (Section 15.31) cannot be accounted for by 
the curvature effect and must be intrinsic. 

In most cases, the X-ray light curve is simply declining early on (possibly with weak 
flaring superimposed), and we observe approximately secularly declining S pea k and HR val- 
ues. In a handful of cases where multiple flares follow a GRB (e.g., Figure El 060124; and 
El 060714), the hardness tracks the flux both upward and downward. Because the brightest 
case (060124 in Figure [3] right) also shows upward and downward Speak trends, we believe 
this behavior is likely responsible for the HR evolution. The parallel or overlapping tracks 
observe d here for bursts with m ultiple flares on the HR-F diagram is also seen for GRB 



pulses ( iBorgonovo fc Ryddl200ll ). 



6. Conclusions 

We have measured the spectral evolution properties for GRBs and afterglows in the Swift 
sample, taken prior to and including GRB 061210. We have established similar spectral 
evolution properties between the X-ray emission coincident with two GRBs (060124 and 
061121) and the X-ray emission in the rapid declines following several GRBs and in 27 flares 
ocurring 10 2 — 10 3 after their GRBs. 

Indirectly from absorbed powerlaw fits which show a time- variable Nh and directly from 
Band model fits, we have derived constraints on the vF v spectrum peak energy S pea k- We 
observe this quantity to evolve in time and to typically cross the XRT bandpass during the 
early X-ray afterglow. Because the X-ray hardness changes little for Band spectra with S pea k 
outsid e the bandpass, the strong hardness variation we observe in >90% of Swift early after- 



glows (IButler fc Kocevskil 120071 ) imply -Ep eak ~ 1 keV, typically. We observe this evolution 
in data taken in both the WT and PC modes (e.g., 050607 and 050714B) and following both 
long duration and short duration (e.g., 050724 and 051227) GRBs. The hardness ratio and 
Speak values scale with the flux as would be expected from the relativistic viewing effects of 
an expanding fireball. This implies that the true variability timescale is even shorter than 
that measured from the observed flare durations. 

Because the late flares are typically softer than the GRB emission, and because the Band 
model a and (3 parameters also appear to evolve in some cases, there is likely an intrinsic 
evolution of the fireball. If the flares are due to shells moving out with lower bulk Lorentz 
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factor or at larger radii than for the prompt emission, we may expect to see differences in 
the ti me properties of flares obse rved at different epochs. This will be explored in a separate 
paper (IKocevski &: Butlerl 120071 ) . If the evolution is occurring on longer timescale at later 
times, when the sensitive XRT is observing, the early X-ray afterglows would provide a 
unique test-bed for theories explaining GRBs, the emission mechanisms, and possibly the 
progenitors. The internal shocks must be active after 10 3 s and must be able to produce 
emission with E peak « 1 keV and very soft (3 <; — 6 (see, also. lZhang et al.ll200q ). Especially 
relevant to the Gamma-ray Large Area Telescope (GLAST), electrons energized by the X-ray 
flares may Compton u pscatter photons at larger radii or in the external shock to the 7-rays 
jKocevski et al.ll2007h . 
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Table 2: 27 Bright XRT Flares 



GRB 


Time Region [s] 


GRB 


Time Region [s] 


050502B 


400-1200 


050712 


150-300 


050730 


130-300 


050730 


300-600 


050730 


600-800 


050822 


410-650 


050904 


350-600 


051117A 


1250-1725 


051117A 


800-1250 


060111A 


200-500 


060124 


300-650 


060124 


650-900 


060204B 


100-270 


060204B 


270-450 


060210 


100-165 


060210 


165-300 


060210 


350-450 


060418 


83-110 


060607A 


93-130 


060607A 


220-400 


060714 


100-125 


060714 


125-160 


060714 


160-230 


060729 


156-300 


060904A 


250-600 


060904A 


600-1000 


060904B 


140-300 







Table 3: -E pea k Evolution Properties 



GRB 


to [a] 


Time Index 


Flux Index 


Data Points Fit 


060124 


510 


-2.9 ± 1.3 


-3.6 ± 1.7 


2 


060124 


555 


-2.2 ±0.3 


-1.8 ±0.5 


3 


060124 


567 


-0.8 ±0.1 


-0.3 ±0.1 


4 


060124 


685 


-2.2 ± 0.4 


-1.2 ±0.2 


3 


060526 


240 


-1.2 ±0.1 


-1.0 ±0.1 


5 


060614 





-2.1 ±0.1 


-0.72 ±0.03 


19 


060729 


75 


-2.0 ±0.5 


-0.4 ±0.1 


5 


060729 


155 


-0.7 ±0.2 


-0.4 ±0.1 


4 


060904B 


140 


-1.3 ±0.3 


-0.7 ±0.2 


5 


060929 


470 


-1.1 ±0.2 


-0.6 ±0.1 


7 



Notes: Changes in the best-fit -E pe ak with time are relative to the start t . The start time is somewhat 
arbitrary, based on the approximate start of each pulse (or flare). 
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Fig. 1. — Selected vF v spectra from GRBs 061121, 060124, and 060614, demonstrating the Band 
model fits to a time varying spectral curvature — as seen in plots with logarithmic axes — and -Speak 
evolution. The X-ray data are corrected for photoelectric absorption using the best fit late-time 
values of Nh in Figures [3] (left). The softest spectrum in the middle panel is divided by a factor 
ten for legibility. The counts spectra are jointly fit by forward folding the Band model through 
the instrument response matrices. For the spectral fits (Table 1), the BAT data are not binned as 
shown here. 
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Time Since BAT Trigger [s] 



Fig. 2.— The light curve for GRB 060124. The X-ray data (0.3-10.0 keV) are plotted in black. 
The shaded regions in the background depict the X-ray light curves in two energy bands (0.3-1.3 
keV and 1.3-10.0 keV) and in the hard X-ray/7-ray bands of BAT (15-100 keV and 100-350 keV). 
The background light curves are each denoised and normalized to their peak intensity. The harder 
regions are darker. The sub-panel shows the early and late XRT light curve. 
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Fig. 3. — Powerlaw (left panels) and Band model (right panels) fits to the GRBs 061121, 060124, 
and 060614. Time-correlated Nh-T variations in the left plots are better modelled by spectral 
models with time-evolving -Break's in the right plots. The Nh values peak when -E pea k ~ 1 keV. 
The powerlaw fits are performed for only the X-ray data, whereas the Band fits (actually nested 
powerlaw then exponential times powerlaw then Band fits, as described in the text) apply to the X- 
ray and 7-ray data. Trends in the Band model parameters, when observed, are fitted and presented 
in Table 3 and in the text. These time variations are given relative t o the a pproximate pulse start 
times. Galactic column densities are taken from lDickey &: Lockmanl (|1990l ) The Band fits use the 
late time Nh values plotted in the left panels, derived from X-ray fits at t > 10 4 s. 
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Fig. 4.— The hard X-ray BAT and XRT light curves for GRB 061121. The late-time light curve 
(a plateau at t ~ 200s followed by a decline beginning at t ~ 3000s) is plotted in the sub-panel. 
The XRT light curve has been multiplied by 5 to bring it above the BAT light curve. 
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Fig. 5. — The hard X-ray BAT and XRT light curves for GRB 060614. The late-time light curve 
(a plateau at t ~ 200s followed by a decline beginning at t ~ 4000s) is plotted in the sub-panel. 
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Fig. 6 — The hardness evolution in GRBs 061121, 060124, and 060614. (top plots) Hardness 
versus rate fit, indexed as a function of time, showing evolution along roughly parallel tracks, 
(middle plots) The X-ray light curve and fit (red curve) as source of the time indexing, (bottom 
plots) The X-ray light curve in each band (hard is red, soft is black) for each time segment and the 
hardness during each time segment. This is well fit during the declines by the square-root of the 
rate (dotted red line) in GRBs 061121 and 060124 and by a power close to the square root of the 
rate for GRB 060614. 
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Fig. 7. — Plots of the hardness and Speak evolution for flares after GRB 060714. The Band 
fits allow only a coarse time resolution, whereas the hardness study demonstrates a fine timescale 
changes in the spectrum which track the flux across flares. -E pea k evolution from Band model fits to 
the BAT and XRT data (top plot). Typical values for the photon indices (a,/3) are given. Hardness 
versus rate fit (middle and bottom plots), indexed as a function of time, showing parallel evolution 
tracks. 
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Fig. 8. — Powerlaw (left panels) and Band model (right panels) fits to the GRBs 060526, 060729, 
060904B, and 060929. Time-correlated Nh~T variations in the left plots are fit by spectral models 
with time-evolving -Speak 's in the right plots. Trends in the Band model parameters, when observed, 
are fitted and presented in Table 3 and in the text. See also Figure El In the Band model plots, a 
values which appear to be above and outside of the plotted range are those which reach and remain 
at the paramater bound a = (see Section I3.1|) . 
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Fig. 9.— Hardness plots for GRBs 060526, 060729, 060904B, and 060929. (top plots) Hardness 
versus rate fits, indexed as a function of time, showing evolution along roughly parallel tracks, 
(middle plots) The X-ray light curve and fit (red curve) as source of the time indexing, (bottom 
plots) The X-ray light curve in each band (hard is red, soft is black) for each time segment and the 
hardness during each time segment. This hardness is well fit during the declines by the rate to a 
power close to 0.5 (dotted red curves). See also Figure 
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Fig. 10. — As also discussed in Paper I, there is an outlier population of very soft Swift XRT 
afterglow time regions with respect to the majority population clustering near photon index T ~ 2. 
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Fig. 11. — Powerlaw fits to high signal-to-noise data (10 4 counts, 0.3-10.0 keV) simulated from 
a Band spectral model with a = — 1 and (5 = —3. Each fit is statistically acceptable (x 2 /^ ~ 1)- 
With the passage of the uF v peak energy E pea ^, the best-fit photon index V steepens smoothly. An 
artificial increase in the inferred X-ray column density Nh linearly proportional to T is observed 
for peak energies -E pea k in the XRT bandpass. The effect is present, with larger error bars, for 
spectra with few counts. 
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Fig. 12. — Time integrated spectral fits to the flares in Table 2, also shown in Figure f!3| demon- 
strate a significant positive correlation between the column density parameter Nh (observed minus 
Galactic) and the photon index T. Although these quantities are correlated for a given spectrum, 
we do not expect a correlation at different times for the same event (see below) or at any time 
for separate events as found here. This is evidence tying the X-ray flares to an excess spectral 
curvature at X-ray wavelengths. 
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Fig. 13. — During the decline phase of the X-ray flares from Table 2 (also Figure fT2j) . the hardness 
ratio (HR), defined as the ratio of counts in the 1.3-10.0 keV band to the counts in the 0.3-1.3 
keV band, correlates strongly (Kendall's tk = 0.6) with the count rate (0.3-10.0 keV), following 
roughly a powerlaw relationship (lef t plot). There is a consistent and long known relation valid for 



a majority of pulses seen in GRBs (jGolenetskii et al.l Il983l ; iKargatis et al 



1995; 



Ford et al 



1995; 



Borgonovo k, Rvde 



2001 



. The hardness also correlates stron gly with the fluence (right plot), as is 



also the case for GRBs ( Liang fc Kargatislll996l ; iRvdd 120051 ) . That is, the hardness evolves more 



rapidly when the flares are brighter. 
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Fig. 14. — Timing statistics for the bright flares in (Table 2). The flare Tqq durations (left plot) and 
rise times (right plot) are systemat i cally l onger in the soft X-ray channel (left plot), by 8± 1% and 
25 ±5%, respectively. iNorris et al.l (|1996l ); lFenimore et al.l (|1995l ); lFenimore. Madras. &: Navakshin 



(1996) discuss similar properties of GRB pulses. 
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Fig. 15. — We believe Njj variations are an incorrect explanation for the spectral evolution in 
the flaring, high— z GR B 050904. These data are coarsely grouped into three time intervals by 
Campana et al.1 (|2006a ) and fit to show a time-decreasing X-ray column density. At finer time 
resolution (left plot), we see that the Njj parameter decreases toward the late time value before 
and after an unphysical increase. The maximum Nh corresponds to -Ep C ak m the XRT band (right 
plot). The hardness during this period tracks the flux to the 0.6 ± 0.2 power (see also, Figure X in 
Butler & Kocevski 2007b), consistent with a Band model spectrum evolving via the curvature-effect 
(Section \5A 
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Table I. Selected Band or Powerlaw*Exponential Model Spectral Fits 



GRB 


Time 

H 


a 


a 


-^peak 

[kcV] 


0.3-10 keV Flux 
[10 -9 erg cm -2 s _1 ] 


xl H 


Signif. 


060124 


569-600 


-1.23 ±0.04 




108l~ 


26.8 ±1.0 


1.21 (154) 


5.9cr 


060124 


400-569 


-1.04 ±0.03 


-2.ol8:S 


27.2±f;| 


8.6i°; 2 


1.01 (410) 


io- 77 


060124 


720-770 


-0 3 +0 - 3 


-2.2 ±0.1 


1.3 ±0.2 


5.8^ 


1.08 (158) 


io- 12 


061121 


60-90 


i i 9 +0.01 
1 - i/: -0.02 




270l£ 


5o.3_ 1 3 


1.07 (270) 


9.3a 


061121 


126-140 


-0 +ao 


-2 4 +ai 

Zl4: -0.2 


0.95+H 5 


3.7 ±0.2 


1.16 (117) 


3.9cr 


060614 


97-111 


-0.7 ±0.1 


-2.4 ±0.1 


8.6 ± 1.2 




0.86 (169) 


1Q -98 


060614 


237-297 


-1.2 ±0.2 


o s+ - 2 
_/ -°-0.3 


1.1 ±0.1 


3.2 ±0.1 


0.97 (246) 


io- 27 



Note. — The quoted errors correspond to the 90% confidence region. The "Signif." column refers 
to the fit improvement significance relative to a simple powerlaw model, determined from a A^ 2 

test. The quoted fluxes are unabsorbed. 



